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STUDIES RELATING TO TEMPERATURE CONTROL 
OF A LARGE SPACE TELESCOPE 

By S. Katzoff 
Langley Research Center 

SUMMARY 

Analytical methods are developed for estimating the circumferential and longitudi- 
nal temperature distributions in a large space telescope, idealized as a simple insulated 
cylindrical tube with a flat mirror across one end. The basic variable is oT^, where 
a is the Stefan-Boltzmann constant and T is the absolute temperature. The circumfer- 
ential distribution is found as a Fourier series. With the view factor between circular 
bands of the tube approximated by a simple exponential, the longitudinal distribution is 
then found as an exponential expression. 

The effects of wall conduction, multilayer insulation, thermal-coating characteris- 
tics, heat pipes, and heated collars are analyzed, with numerical examples. The effects 
of nonuniform coatings and insulation are also discussed. A brief analysis is also given 
of unsteady heat flow through the insulation, as in alternate light and shadow of a low orbit. 
The numerical examples were mostly simplified by assuming that the only external ther- 
mal input was from solar radiation directed at right angles to the tube axis. 

Multilayer insulation around the large tube can limit the circumferential tempera- 
ture nonuniformity to tenths of a degree. Adding circumferential heat pipes can reduce 
it further, to hundredths of a degree. The insulation, however, lowers the tube tempera- 
ture and causes a pronounced drop in temperature along the tube toward the open end. A 
heated collar around the inside of the tube near the open end can overcome this problem 
and provide a nearly uniform temperature along the tube between the mirror and the 
collar. 


INTRODUCTION 

Among the more ambitious of the space -science proposals that are under study by 
various research groups is the diffraction -limited Large Orbiting Telescope (LOT), or 
Large Space Telescope (LST), more or less as it was proposed by Lyman Spitzer, Jr., 
some years ago. (See ref. 1.) The extent and intensity of interest in such a telescope 
are attested by the many papers dealing with the various facets of its design and con- 
struction that were presented in a 1969 workshop on optical telescope technology. (See 
ref. 2.) 



As was pointed out in reference 1, maintaining diffraction -limited optics in a tele- 
scope of, say, a 3 -meter or larger diameter, implies an extremely fine temperature 
control in order to avoid thermal distortions. Temperature uniformity and constancy to 
within a small fraction of a degree may be required. However, the availability of mate- 
rials with very low thermal coefficients of expansion should help to alleviate the problem 
Active optics (that is, continuously adjustable mirror figure and alinement) is an addi- 
tional approach to alleviating the problem. (See ref. 3.) Heat pipes, vapor chambers, 
phase -change materials, thermal shields and enclosures, multilayer insulations, and 
automatic shutters may also find application in solving the thermal problem. 


The present paper is an analytical study of the thermal problem for a simple ideal- 
ized configuration that should serve as a first approximation to any eventual design. Spe- 
cifically, the configuration is essentially a cylindrical 
tube with its inner surface perfectly black, and with a 
flat, perfectly reflecting mirror at its closed end. (See 
sketch (a) .) The axis is assumed to be normal to the 
solar radiation for most of the analyses. Although 


•Mirror 


♦ I 


♦ ♦MM 

Sketch (a) 

the large telescope would probably be of the Cassegrain type, no effort has been made to 
take into account the presence of the secondary mirror or of the hole at the center of the 

primary mirror. This simplification should not greatly affect the usefulness of the anal- 
ysis and results. 


This study has two purposes: (1) to present analytical methods that are available 
for making such a study, and (2) to present enough numerical results to indicate the 
degree of temperature uniformity attainable by the various means discussed. The ana- 
lytical methods concern the determination of circumferential and longitudinal temperature 
distributions in the tube wall. The analysis of the circumferential distribution involves a 
Fourier -series approach; the analysis of the longitudinal distribution involves use of an • 
exponential view factor between circular bands in the cylinder. None of the methods 

described require more computer capability than is available in a programable desk-top 
calculator. 

Steady-state conditions were assumed for nearly all cases. As previously indicated, 
for most of the calculations the cylinder axis was taken as normal to the solar rays, and 
the shaded half of the cylinder was assumed to receive no radiation at all. The effects of 
insulation, thermal coating, heat pipes, and heated collars were considered. A few cal- 
culations concerning the effect of alternating sunlight and shadow, as in a typical low orbit, 
were also made. No effort was made to cover wide ranges of the parameters and of their 
various combinations. However, the results should be adequate for showing the main 
aspects of the problem. 
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Some of the results given in this paper were presented previously in a paper at the 
workshop (ref. 4). The analytical methods, however, were not described in that paper. 

This paper should be recognized as fairly limited in both purpose and scope. It 
mainly concerns the temperature of the telescope tube. If this temperature can be main- 
tained uniform and constant, the uniformity and constancy will be reflected in (1) the 
radiation striking the face of the primary mirror, (2) the radiation environment of the 
secondary mirror, and (3) the alinement of the secondary mirror with the primary mirror. 
But even within this limited purpose of the paper, some important pertinent problems have 
not be analyzed - for example, the use of a removable cover for the tube opening. 

SYMBOLS 

constants in expressions for z(x) and y(x) (for finite -length tube) 

circumferential mean value of z on outer layer of insulation (for finite - 
length tube) 

coefficient of cos n0 in Fourier series for z 
mean value of z around infinite cylinder 

solar constant, ~ cal/cm 2 -sec , 

thermal conductivity of wall material, cal/cm-sec-K 


function defined by equation (6) 


H + Pi c 

h = pl - 

VI + g 

k "radiation conductance" of multilayer insulation (Rate of heat transmission 

is k(z Q - z A ), cal/cm 2 -sec.) 


A,B,C 


b 

C 


c = 


ke 

k + e 


cos* 6 




Ce^ 
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q = oT 4 
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T 

T 
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u 

X 

X 

y = oT 4 

4 

z = ctT 

a 

6 

e' 


length of telescope tube (mirror to opening), radii 
distance from mirror to beginning of heated isothermal collar, radii 
distance from mirror to end of heated isothermal collar, radii 
half the spacing between circumferential heat pipes, cm 
coefficients in power series for z(x) 
on heated isothermal collar 
radius of cylinder, cm 

factor for radiation exchange on inner wall of finite -length tube 
temperature, K 

circumferential root mean fourth power of wall temperature, ^ ^ , K 
time, sec 

rate at which radiation leaves a unit area, cal/cm 2 -sec 
longitudinal distance along tube from mirror, radii 

longitudinal distance from circle midway between adjacent heat pipes, cm 
along tube between heated isothermal collar and open end 

solar absorptivity of outside surface 
angular distance around cylinder from subsolar point 
dummy variable corresponding to 6 
dummy variable corresponding to x 
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Pi 


reflectivity of inner cylinder wall 


e thermal emissivity (of outside surface unless otherwise indicated) 

a Stefan -Boltzmann constant, 1.3556 x 10'^ cal/cm^-sec-K^ 

t wall thickness, cm 

Subscripts: 

0 outside surface of insulation 

1 tube wall 

f finite -length tube 

m midway between circumferential heat pipes 

max maximum around circle 

min minimum around circle 

1 outside surface of variant section (see sketch (j)) 

0,1, 2,. ..,9 ten sheets of a multilayer insulation (eq. (34)); 9 refers to sheet adjacent 
to tube wall, 0 refers to outside sheet 

10 tube wall 

Primes denote differentiation with respect to x. 

INFINITELY LONG CYLINDER 

As a first step, the case of the infinitely long thin -walled metal cylinder will be 
analyzed. The cylinder is in space, at one astronomical unit (1 A.U.) from the sun, with 
its axis normal to the solar radiation. No other radiation is present, so only the sunlit 
half receives radiation, while the other half receives no radiation at all. 

The following section will consider the simple case of the uninsulated cylinder in 
order to exemplify the method of analysis. Subsequent sections will consider the effects 
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of (1) insulation, (2) circumferential heat conduction in the cylinder wall, and (3) circum- 
ferential heat pipes attached to the cylinder wall. The analysis is similar in many 
respects to that of reference 5. 


Uninsulated Cylinder 

Internal radiation exchange .- The inner surface of 
the cylinder is assumed to be perfectly black. Then, of 
the radiation emitted from any (infinitely long) longitu- 
dinal element of the internal surface at 9 1 , the fraction 

i le' el 

- sin - — s — - de is received by an element at 6 of 
4 ^ 

angular width d 6 (ref. 6, pp. 19-21). (See sketch (b).) 

The basic variable will be taken, not as tempera- 

4 

ture T, but as z = oT , where a is the Stefan- 
Boltzmann constant. Let the distribution of z around 
the cylinder be given by a cosine series in 9, the angle 
from the subsolar point: 

z(0) = oT 4 (0) = a Q + a 1 cos 9 + a 2 cos 29 + a g cos 3 9 + . . . (1) 

Then at any circumferential location 9, the received radiation intensity contributed by 
the cos n0 term is 

i n8+2ir r' , a n 

7 a \ cos n0 sin ■■ ■■ 0 d 8 = — cos n0 (2) 

4 n «V=0 2 4n2 - 1 

Thus, the received radiation intensity also varies as cos n 8, but its phase is reversed 
and its amplitude is reduced by the factor 4n^ - 1. 

Temperature distribution .- Let the cylinder diameter be so large relative to the 
wall thickness (for example, a 3 -meter -diameter cylinder with a wall of 1 -mm -thick 
aluminum) that heat transfer by thermal conduction around the wall is negligible in com- 
parison with heat transfer by radiation across the inside of the cylinder. 

Let 



b solar constant, cal/cm2-sec 

Oi solar absorptivity of external surface 

e thermal emissivity of external surface 
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As already noted, the internal surface is assumed to be a perfect absorber and emitter. 
Then, on the sunlit half of the cylinder^) 0 |< 


ab cos 


GO 00 

0 - y — — cos n0 = (1 + e) ) a n 

. A 4n 2 - 1 A 


cos n0 


I) <3) 


where the first term on the left is the heat input from the sun; the second term on the left 
is the heat received by radiation from the entire internal surface of the cylinder (see 

eq. (2)); the term e ^ a n cos n0 on the right is the heat radiated away from the outer 

surface; and the remaining term on the right, ^ a n cos n 0, is the heat radiated from the 
inner surface. 

For the shaded half of the cylinder < | 0 | < 77 ^, the first term on the left-hand side 
of equation (3) is omitted, and the equation becomes 


OO OO 

_\ — cos n 0 = (1 + e) ) a n cos n 9 

^ 4n 2 - 1 V 




Equations (3) and (4) may be combined by 
introducing the function cos* 9, which is the 

same as cos 6 for | 6 |< and which is zero 

for - < | 9 ) < 77 . (See sketch (c).) The combined 



equation is 


0 = 0 ?r/2 77 


Sketch (c) 


OO 00 

ab cos* 9 - ) — — cos n6 = (1 + e) ) a n cos n0 
V 4n 2 - 1 ^ 


0 < 0 < tt) ° (5) 


The function cos* 0 may be expanded in a Fourier series: 

cos* 0 = - + + JL cos 20 - rr— cos 40 + cos 60 - 

tt 2 3-ir 157 t 3577 


, / i\2 2 cos n0 a 

+ 1-1) 7-5 r- 0 

(n 2 - 1)77 

Substituting this series for cos* 0 in equation (5) gives 


(n even) (6) 


. 1 cos 0 2 

ab - + — - — + ■=— cos 20 - 
2 377 


+ 1 + e ]a n cos n0 
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The a n are now found from this equation simply by equating coefficients of cos n0. 
The results are 



With these coefficients, the series (1) may be evaluated to give oT 4 (0), from which 
T(0) may be determined. The constant term of the series, ao, which is the circumferen- 
tial average value of oT 4 , depends only on the a/e ratio, and not on the separate values 
of a and e . 

Figure 1 shows the calculated circumferential temperature distributions for the fol- 
lowing four pairs of values of external absorptivity and emissivity: 


Case 

a 

€ 

1 

1.0 

1.0 

2 

.4 

.4 

3 

.25 

.85 

4 

.15 

.85 
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The third and fourth cases correspond roughly to white paints. (One might also consider 
the fourth to represent a fresh white paint, and the third to represent the same paint after 
extended exposure in space.) The second case might correspond roughly to a metallic 
paint, like fine aluminum flake in a suitable clear vehicle. The first case, of course, 
corresponds to a theoretically black surface. 

The root mean fourth power of the temperature (ag/a) 1 / 4 , designated as T, has 
been indicated by a short line for each case. Its value is somewhat higher than the true 
average temperature. 

The curves have about the expected shapes, but several features might be mentioned. 


(1) The coolest part of the wall is not at the antisolar point, 0 = 180°, but at 
6 = 90°. The reason is that the hot subsolar area around 6 = 0° radiates more strongly 
to the antisolar point than to any other point, because of the sin — — factor, which is 
a maximum when | 0 f - 0| = 180°; whereas a point just at the edge of the shadow receives 
the least radiation. As will be shown subsequently, circumferential heat conduction in 
the wall will eliminate the sharp breaks in the curves at 0 = and will also shift the 
location of the temperature minimum to larger angles. 


(2) The curves all seem to pass through T near 0 = 70°. 


(3) Although cases (1) and (2) have the same value of a/e , the first shows a much 
larger temperature variation than the second. The reason is that the temperature - 
equalizing tendency of the internal radiation exchange is the same in both cases (since 
the internal surface was assumed to be black) , and is thus relatively more important in 

the second case, where the external heat exchange is weaker (that is, lower a and e) 

than in the first case. 


(4) The largest temperature difference between the warmest and coolest points is 
99 K (for a = e = 1.0); and the smallest temperature difference is 58 K (for a = 0.15; 
e = 0.85). Differences of this order are unacceptable for the telescope. It is apparent 
that a very effective insulation will be essential in order for the rate of heat flow between 
the outside and the inside to be very small relative to the rate of heat exchange by radia- 
tion within the cylinder. 

The e ffect of a gray internal surface . - In the preceding analysis the internal sur- 
face was assumed to be perfectly black, for which case the fraction of the radiation 
emitted from an element at 0' that is received by an element of angular width d0 at 
location 0 is ^ sin 1 ® ~ d 0. If the internal emissivity is some value q < 1, not 
only is the emitted radiation reduced by the factor ei (forgiven T(0')), but the fraction 
of this radiation that is absorbed by the element at 0 (including what is absorbed after 
one or more reflections) is changed to 
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de 

(See ref. 7; a related analysis is given in ref. 8.) For values of about 0.90 to 0.95, 
the value of this expression is less than i sin ! -£■ ' - 6 1 for le' - e! > 81°, and exceeds 
it for 1 0’ - fl|< 81°. 2 

Correspondingly, equation (2) for the absorbed radiation intensity contributed by 
the cos n$ term is changed to 




If the factor 


— is substituted for 

4 „ 2 - «i 


cients (eqs. (8)) become 


— in equation (7), the Fourier coeffi- 

4n 2 - 1 



where the outside emissivity has been written e 0 to distinguish it from the internal 
emissivity e^. 

Values of e ^ of about 0.95 have been reported for some mat black coatings such 
as would presumably be used in the telescope tube. With such values, the temperature 
variation around the tube would be increased (relative to what is shown in fig. 1) by 
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about 5 percent. Accordingly, the results shown in the present paper, which are all 
based on ej = 1, slightly underestimate the circumferential temperature variation. 


The Effect of Insulation 

Assumed insulation characteristics.- Consider the cylinder to be wrapped in multi- 
layer insulation — that is, many layers of aluminized plastic sheet, loosely packed. The 
multilayer insulation is idealized in the sense that heat transfer across it will be assumed 
to occur only by radiation between adjacent layers. Such idealization provides important 
simplification of the analysis because it permits z(soT 4 ) to be retained as the basic 
variable in all parts of the analysis; in particular, the heat flow through the insulation is 
proportional to the difference between the z values on the outer and inner surfaces. 

Thus, the heat flow is written as k(z Q - Zj), where k is a "radiation conductance,” z 0 
is the value of z on the outermost layer of insulation, and zi is the value of z on the 
inner cylinder. For a small difference between the temperatures T 0 and T^, the heat 
flow would then be 

kA(oT 4 )«4koT 3 AT 

where the factor 4koT 3 may be interpreted as the usual conductance, which, however, 
varies rapidly with T. 

In order to test the approximate validity of this concept, one might examine whether 
measured conductivities of multilayer insulations do indeed vary with T 3 . In refer- 
ences 9 and 10 (p. 10) are a set of data for several different multilayer insulations, obtained 
with a very sensitive apparatus that could provide conductivity values even when AT 
was as small as 6 K. These data have been replotted on a log-log plot in figure 2. Best- 
fit straight lines are shown for the various insulations; in addition, an exact T 3 varia- 
tion (the heavy solid line) has been added to the plot. It will be seen, first, that most of 
the groups of data points do indeed lie nearly on straight lines; and furthermore it will be 
seen that the slopes of most of the lines lie between 2 and 3, whereas two of the lines 
have slopes exceeding 3 (3.5 and 4.25). The high slopes are not necessarily in error, 
since theory predicts that the aluminum coatings should have an emissivity approximately 
proportional to T (as may be shown, for example, from eq. (2) of ref. 11, by using hand- 
book values of the electrical resistivity of pure aluminum), which would cause the con- 
ductivities to vary as T 4 rather than as T 3 . The lines with the smallest slopes are 
for the insulations with the highest packing densities, for which there is presumably an 
increased proportion of conductive heat transfer. On the whole, then, one may conclude 
from figure 2 that the T 3 variation assumed for the present analyses is not unreason- 
able for a well constructed multilayer insulation that is not too densely packed. 

In the calculations to be presented later, k values of 0.001, 0.0025, and 0.025 
were used for the insulation. Since the thermal emissivity of aluminized Mylar sheet is 
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about 0.05, these values correspond theoretically to about 25, 10, and 1 or 2 layers, 
respectively. The heavy solid line of figure 2 corresponds to a 1 -cm -thick insulation 
for which k = 0.0025. Since the line lies within the range of test results, values of the 
order of 0.0025 may thus be considered as fairly realistic (although it is well known that 
careful design and construction are necessary to prevent extraneous heat flows that 
greatly exceed predictions from such laboratory values). 


Derivation of Fourier coefficients . - The problem is set up as before, although it is 
now slightly more complicated. On the sunlit half of the cylinder, the heat flow through 
the insulation satisfies the equation 

o-b cos 0 = ez 0 + k(z 0 - Zi) (l^f) ( 9 ) 


where the left-hand side of the equation is the solar heat input to the outermost sheet of 
insulation; the first term on the right is the heat radiated to space from this outermost 
sheet; and the second term on the right is the heat transferred through the insulation 
from the outermost sheet to the inner cylinder. The z terms are functions of 8. On 
the shaded half of the cylinder, the solar heat input is zero, and the equation is 

0 = ez 0 + k(z Q - Zi ) (| < I 0 I < 71-) (10) 

Equations (9) and (10) may be combined, as before, by writing 

ah cos* 8 = ez Q + k(z Q - zj\ (11) 

Solving for z Q gives 

ab cos* 8 + kzj 


whence 


" z i 


ah cos* 8 - ezi 
e + k 


( 12 ) 


Equilibrium of the inner cylinder provides the equation 





cos nd = 


l 


l n 


cos n0 


(13) 


where the first term on the left is the heat transferred in through the insulation at 8, 
the second term on the left is the radiation received at 8 from the entire inner surface 
of the cylinder, and the term on the right is the heat radiated away from the element at 
8 into the inside of the cylinder. Substituting the expression for z 0 - z j from equa- 
tion (12) into equation (13), while replacing cos* 8 with its Fourier series and replacing 

oo 

z i with ^ a n cos n0, gives 
0 
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kg?b /. 1 cos 0 + 

6 + k \7T 2 


~ cos 20 - 
37T 



a n cos n0 


from which, by equating coefficients of cos n0, the Fourier coefficients for z^ are 
obtained: 


a, = 


ab 




2ab 

a2 " 3 TrfM ejjc 

15 k 


(14) 


The coefficient a Q (the average value of zi) is the same as before. The other 
coefficients are changed only by the factor (e + k)/k on one of the terms in the denom- 
inator. With k = °°, the coefficients reduce to those previously derived for the uninsu- 
lated cylinder (eq. (8)). 


The Effect of Circumferential Heat Conduction 


It was previously indicated that for the Large Space Telescope, the heat conduction 
around the cylinder would be negligible relative to the radiation heat transfer. It would 
not be negligible for a much smaller telescope, however; so, for completeness, the pre- 
ceding analysis will be extended to include this effect. The equations will be developed 
for the general case in which the wall thickness is r cm and the conductivity of the wall 
material is C cal/cm-sec-K. In all the computed examples, however, the wall will be 
assumed to be 0.1-cm -thick aluminum (t = 0.1; C = 0.5; Ct = 0.05). 

The circumferential flow of heat, by conduction, into 1 cm 2 of the wall is 

4 *^ z i 3 dT 

where r is the cylinder radius. With za = oT , — — = 4oT — , or, since T 

_ d0 d0 

close to its root mean fourth power T, 


Ct 


d 2 T 

r 2 d0 2 ’ 


remains 


~ 4oT 3 — 
d0 d0 


d2z i -3 d 2 T 
1 —s 4oT' j a 1 


d0^ 


d02 
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The Effect of Circumferential Heat Pipes 

Derivation of heat -transfer term .- Suppose that the infinite cylinder has circumfer- 
ential heat pipes around it at small fixed intervals. (See sketch (d).) These pipes would 
serve to reduce the circumferential temperature variation. For 
the present analysis it will be assumed that the heat pipes are 
perfect conductors and that the technique of bonding them to the 
wall provides negligible thermal resistance between the pipes 
and the wall. They thus provide a periodic distribution of iso- 
thermal rings along the cylinder. Admittedly, present tech- 
nology frequently produces only a poor approximation to the Sketch (d) 
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assumed ideal conditions. However, it may be hoped that the present research and 
development effort in this field will result in significant advances in the art during the 
next few years. 

The heat pipes acquire the average temperature of the cylinder. For most of the 
sunlit half, heat is delivered to the heat pipes, flowing along the wall toward them; and in 
the shaded half, heat is removed from the heat pipes, flowing along the wall away from 
them. Correspondingly, along any longitudinal cylinder element there will be a periodic 
variation of temperature, about as shown in sketch (e). The short arcs in the sketch are 

Heat-pipe locations J ghaded area 



Sketch (e) 

approximately parabolic. (In the one -dimensional case with uniform addition or loss of 


= Constant, for which the 


d^T 

heat by radiation, the governing differential equation is — *- 

cbsr 

variation of T with x would be exactly parabolic.) 

Two modifications to the analysis result: 

(1) Corresponding to the longitudinal heat flow in the wall, the rate of loss of heat 

2 

from a cm 2 of wall is -Ct -^ j-, where, as before, the factor Cr is the product of wall 

ax 2 

thickness (r cm) and conductivity (C cal/cm-sec-K). The parabolic variation gives 

s2l 


AT = (AT) 


m 


‘-'I 


from which 


where 


9 2 T 

ax 2 


2 (AT) 


m 


P 

AT 

AT m 


half the spacing between heat pipes 


wall temperature minus heat -pipe temperature 


value of AT midway between adjacent heat pipes 


longitudinal distance from midway point, cm 
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Since the range of variation of T is small, the parabolic distribution of T implies 

00 

a parabolic distribution of z (= oT 4 ) also. Let z i = a 0 + ^ a n cos n0 be oT 4 on the 

1 

circle that is midway between adjacent heat pipes, where a Q , the average value of z i; 

oo 

is also the heat pipe value of oT 4 . Then AT m corresponds to ^ a n cos n 6; or, by the 
relationship derived in the preceding section, 4 

Az v 


AT m = 


J m 




cos n6 


4oT 3 4a 1//4 a 0 3//4 \ 

The conduction heat flow out of the wall toward the heat pipes, per cm 2 , is then 


-Ct^T= -Ct( 


-2 AT 


ml 


Ct 


OO 


2 ~ l/4„ 3/4 Z a n C0S ne 


2p*o L ' *a 0 -- i 


(2) Since the average height of a parabolic arc is two -thirds of its maximum height 

OO 

(sketch (f)), the average z along any infinite longitudinal element is a Q + |^ a n cos n 6. 
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This average value will be used in the term for the radiation received 
at 6 from the inside of the cylinder (a modification of eq. (2)). 



Sketch (f) 


It must be noted that the assumed fairly small heat pipe spacing is necessary for 
the preceding discussion to be valid. For example, if the spacing is very large, the lon- 
gitudinal temperature distribution would no longer be approximately parabolic, so that 
neither the two-thirds factor nor the expression for 9 2 T/ax 2 would be valid. 


Heat -balance equation and solution . - The heat -balance equation, including the terms 
for insulation, circumferential conduction, and longitudinal conduction to heat pipes, is, 
for the circle that is midway between adjacent heat pipes, 




— = a n cos n0 

4n 2 - 1 
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( 16 ) 


where the first term on the left is the heat transmitted in through the insulation; the 
second and third terms on the left represent the radiation received from the inside of the 
cylinder; the first two terms on the right represent the heat radiated away into the inside 
of the cylinder; the third term on the right is the heat conducted away circumferentially; 
and the fourth term on the right is the heat conducted away longitudinally. 
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As before, the factor (z Q - z^) is replaced by the expression in equation (12); and 
the Fourier coefficients are then determined by the procedure previously described. 
They are 
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These coefficients reduce to those for the case without heat pipes (coefficients (15)) if the 
last term in the brackets is deleted and if the factor 2/3 is changed to 1. If the fourth 
term in the brackets is also deleted, the coefficients (14) are obtained. Finally, if k is 
made infinite, the coefficients (8) are obtained. 


Calculated Results 

A number of circumferential temperature distributions were calculated by means 
of the sets of Fourier coefficients (14), (15), and (17). As previously mentioned, for all 
cases in which heat conduction in the wall was considered, the product of wall thickness 
and conductivity was taken as 0.05. 

The cylinder radius r was taken as 150 cm for most of the calculations, but a few 
calculations were made for a radius of 15 cm. The effect of circumferential conduction 
(given by the term containing r in the coefficients (15) and (17)) is not important for 
r = 150 cm; but since r occurs as the inverse square, the effect rapidly increases in 
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importance as r decreases. The heat -pipe spacing 2p was taken as 40 cm for all 
heat -pipe calculations. The isothermalizing effect of heat pipes also varies inversely 

O 

as the square, as indicated by the p factor in the last term in the brackets in the 
coefficients (17). 

Figure 3 shows circumferential temperature distributions calculated with circum- 
ferential conduction taken into account. The cylinder radius was taken as 150 cm for all 
cases but one. For the insulation, k values of 0.001, 0.0025, 0.025, and °° were used 
(k = °° corresponds to no insulation). In figure 3(a), for k = °°, a curve from the zero- 
conduction results of figure 1 has been added (the dashed curve). It is seen that the cir- 
cumferential conduction has relatively little effect for a large cylinder (r = 150 cm), 
although it does serve to wash out the break at 8 = 90°. On the other hand, when r is 
only 15 cm, the temperature -equalizing effect of conduction is seen to be very pronounced, 
since the circumferential variation is now only half as much as when r = 150 cm. 

Comparison of the curves in figure 3(a) with the upper set in figure 3(b) shows that 
even a little insulation will greatly increase the temperature uniformity. For k = 0.025, 
corresponding to only one or two layers of aluminized plastic sheet, the variation was 
reduced by a factor of about 15. For k = 0.0025, the variation is only a few tenths of a 
degree. 

Comparison of the three sets of curves in figure 3(b) shows that for a given a 
and e, the temperature variations are approximately proportional to k. In fact, for 
the two smallest values of k (the lower two sets of curves), the proportionality to k 
is almost exact. 

The curves in figure 4 show the further improvement in temperature uniformity 
that can be obtained with heat pipes. The upper set of curves, for the zero -insulation 
case, shows improvements by factors of 1.9 to 3.3, compared with the curves of fig- 
ure 3(a). The lower set, for k = 0.0025, shows improvements by factors of 2.3 to 5.7. 

As already indicated, reducing the heat -pipe spacing would cause rapid further improve- 
ment. It must also be noted that these results in figure 4 are for the circle midway 
between adjacent heat pipes. The average variation around the cylinder is only about 
two-thirds as much. 

Calculations were also made for k = 0.025 and 0.001 for the heat -pipe case. The 
results are not shown, since the curves are again simply proportional to k, and may be 
derived from the results shown in figure 4 for k = 0.0025. 

The quantities of heat transported along the heat pipes turn out to be very small. 

For example, for the middle curve of figure 4(b) (a- = 0.25, e, = 0.85, k = 0.0025, 

AT m = 0.052 K), the heat delivered to each heat pipe between 6 = 0° and 0 = 70° (and 
withdrawn from each heat pipe between 8 = 70° and 8 = 180°) is only 0.061 cal/sec. 
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Accordingly, the original assumption that the heat pipes are essentially isothermal seems 
to be justified. The other assumption — that the bonding between the heat pipes and the 
wall has negligible thermal resistance — seems, with present technology, to be more 
questionable. 


FINITE -LENGTH CYLINDER 

The preceding analyses of the infinitely long cylinder show how to calculate the 
degree of circumferential temperature nonuniformity for given wall, insulation, thermal 
coating, and heat -pipe spacing. The present section will consider the longitudinal tem- 
perature variation along a cylinder of finite length. It will be first assumed that, by 
means of insulation and heat pipes, the circumferential temperature variation has been 
made negligibly small or, at least, has been made so small that it can be considered 
independently of the longitudinal temperature variation. That is, the internal wall tem- 
perature will be considered as a function of longitudinal location only. Some discussion 
and analyses of circumferential temperature variation will be given later. 

Since circumferential heat flow is not considered for the present, it will be con- 
venient, in the analysis, to take the cylinder radius as unity. 

Net Heat Flow Through the Insulation 

In the preceding analyses the average value ag of oT 4 , either outside or inside 
the infinite cylinder, was ab/ire. In a finite -length cylinder, the internal oT 4 is 
everywhere less than this value because heat is lost by radiation out of the open end. 
Correspondingly, the average external oT 4 at any longitudinal location is less than 
ab/ire, because the outside is the source of this lost heat and is continuously supplying 
it through the insulation. The symbol a 0 will be retained, for the present, for ab/ire, 
and the symbol a will be used to represent the local circumferential average oT 4 of 
the outermost layer of insulation. The local circumferential average oT 4 of the cyl- 
inder wall will be designated z. 

It is apparent that the local rate of heat supply through the insulation is propor- 
tional to a - z. However, it is also proportional to a 0 - z. The proof follows: At any 
given longitudinal location along the (unit -radius) tube, the heat -flow equation for a cir- 
cular band 1 cm wide is 

2ab = 27rae + 27r(a - z)k 

where the left-hand side represents the rate of absorption of solar heat by the outer layer 
of insulation; the first term on the right is the rate of radiation of heat away from the ■ 
outer layer; and the last term is the rate of transfer of heat through the insulation to the 
inner cylinder wall. Solving for a gives 
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_ ah + 7rkz 
7r(k + e) 

Substituting this expression for a in the term 27rk(a - z), which is the rate of heat 
transfer in through the insulation, gives 


27rk(a - z) = 


27rk 


«b + 7rkz 
7r(k + e) 



= 27rkj 


ob - Tiez 

7r(k + e) 


27rke fab _ A 
k + e We 1 


= 27rc^aQ - z'j 


(18) 


kc 

where c is . Thus, the rate of heat transfer in through the insulation is propor- 

k + e 

tional to ag - z, as was to be proved. The product c(ag - z) is the average heat trans- 
ferred in per unit area at that location, and will be used in this form in the subsequent 
analyses. 


The Exponential View Factor 


Given a unit-radius cylinder (sketch (g)) with its inner surface perfectly black, and 

given two elemental circular bands of this surface, ■■ ■ 

one at the origin, the other at x, the fraction of the 
total radiation emitted from the first band that is — 
transmitted to the second band, per unit width of r = 1 

the second band, is (ref. 12) -1* 


1 

2 


1 - 


|xl(x 2 + 6)1 



Sketch (g) 


It was shown in reference 12 that this view factor is roughly approximated by the simple 

1 - lx! 

exponential expression — e ■; and that use of the exponential expression greatly simpli- 

Li 

fies the analysis of radiation exchange within a cylinder. 


The two functions are compared in figure 5. The areas under the two curves are 
the same. The expressions have the same values at |x| = 0, 1.106, 4.608, and «>. At 
large |x|, the exponential curve lies below that for the perfect radiator. One might hence 
argue that the exponential curve is the more representative of actual mat black finishes, 
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which tend to become reflectors for radiation received at glancing angles. On the other 
hand, the exponential curve lies above the curve for the perfect radiator for x between 
1.106 and 4.608, and hence, in this respect, it represents a thermodynamically impossible 
surface finish. Of course, if the exponential curve is uniformly lowered by an emissivity 
factor of, say, 0.85, it would lie nearly completely below the curve for the perfect radia- 
tor, so that this objection would not apply. However, 0.85 is considerably below the emis- 
sivities of available black finishes. A better approach might be to use a sum of two or 

three exponentials (for example, 0.605e~*"^ i * * * * * * * * x l - 0.105e~3’^l x ^ to provide a better 
approximation to the curve for the perfect radiator, as was suggested in references 12 
to 15. The mathematics remains tractable, although, of course, it becomes much more 
cumbersome. However, the slightly improved analytical accuracy attainable with such 
modifications would have little practical significance for the present study and would not 
justify the greatly increased computational effort. Accordingly, the simple exponential 

view factor - e 1 1 will be used in the present analyses. Some additional analytical 

justification, other than the similarity of the two curves of figure 5, will be given later. 

A practical consideration is simply that the paint and insulation characteristics would not 
be so accurately known or so accurately uniform and reproducible that extreme effort to 
achieve analytical accuracy is warranted. Furthermore, since the assumed configuration 
does not take into account the presence of the secondary mirror and its supporting struc- 
ture, the curvature of the primary mirror, or the hole in its center, it would hardly be 
reasonable to insist on ultimate accuracy in the analysis. 

The Simple Open Tube 

Consider an insulated, open-ended, internally black tube of unit radius and length 
21, with its axis normal to the solar rays. (See sketch (h).) For present purposes, it 

i 

Mirror 

I 

x = -l x = 0 x = l 

Sketch (h) 

represents an idealized telescope tube of length l with a perfect flat mirror at one end 

(x = 0). The mirror reflects the inside of the tube and thus, in effect, provides for the 

analysis a tube of length 21 with both ends open. 

Use of the exponential view factor provides the following integral equation for the 

heat balance of a unit area at any longitudinal location x: 
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( 19 ) 


z(x) = c[a 0 - z(x)] + ± J* z(|)e^ -x d^ + | ^ zU)e x ^d£ 

— Z x 

where | is a dummy variable corresponding to x. 

The term on the left represents the heat radiated away from the inner surface at x, 
per unit area; the first term on the right represents the heat received at x, per unit area, 
from the outside, through the insulation (see eq. (18)); the first integral represents the 
radiation received from the part of the wall that lies to the left of x; and the second 
integral represents the radiation received from the part of the wall that lies to the right 
of x. 

Transfer the first term on the right to the left-hand side of the equation to get 

(1 + c)z(x) - a Q c = i j . ' . . + i J . * * (19a) 

Differentiate with respect to x to get 

(1 + c)z'(x) z(0 e^" x d^ -|>jstx) +| J z(t) e x "%4 

“Z x 

Differentiate again to obtain 

(1 + c)z”(x) = - |z(x) +| J z(S) e^" x d£ - |z(x) +| J z(£) e x "^d£ 

— Z x 

The two integrals on the right are now the same as the original ones in equation (19a); 
therefore, they may be replaced by the left-hand side of that equation. Thus, a differ- 
ential equation is obtained: 

(1 + c)z" = -z + (1 + c)z - caQ 

which simplifies to 

(1 + c)z" = c(z - ag) 

the solution of which is 

z = aQ - A cosh fx (20) 

where f = \jc/(l + c) and the coefficient A is to be determined. Since z is obviously 
symmetrical with respect to the origin, there will be no sinh fx term in the solution. 
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In order to determine A, substitute z from equation (20) into equation (19a) and 
perform the integrations for some convenient value of x — say, for x = 0.^- Thus, 

(1 + c)(a 0 - A) - a Q c = | J ( a 0 ' A cosh fx ) e ^ d £ + f J Q ( a 0 “ A cosh f ^) e ^ 

where the two integrals are obviously equal. The integration proceeds without difficulty 
and provides a simple linear equation for A, the solution of which is 


A = 


2 a, 


0 


a~U 


Jl 


( 21 ) 


1 +f 1 - f 


Putting this expression for A into equation (20) finally gives the equation for z as a 
function of the distance in radii along the cylinder x. Since z = oT^, the temperature T 
along the cylinder is then easily obtained. 

Calculations were made of the temperature distribution along the cylinder for the 
following listed values of the basic parameters: 

Z: 6, 8, 12, 20 radii 

a: 0.15] 0.25 

e: 0.85), 0.85 

k: 0.001, 0.0025, 0.025, °° 

In all cases the cylinder axis was assumed to be normal to the solar rays, the shaded half 
receiving no external radiation input at all. 

The calculated results are shown in figure 6. It will be seen that for the uninsulated 
cylinders (k = *>), the temperatures remain nearly equal to those for the infinitely long 
cylinders, except near the open end, where they drop sharply. For the insulated cylinders, 
the temperatures even at x = 0 (that is, at the mirror) are well below those of the cor- 
responding infinite cylinders, and they drop steadily toward the open end. 



The Effect of a Gray Internal Surface 

Any real surface has. an emissivity less than 1. The radiation from a surface is 
thus everywhere less than the local z. Furthermore, where this emitted radiation 
encounters some other part of the wall, it is only partly absorbed, and the remainder is 
reflected; and where the reflected part meets the wall it is again only partly absorbed, 

1 An alternative approach is to perform the integrations without giving x a par- 
ticular value, and then equate coefficients. In more complicated problems, however, this 
approach may lead to much difficult algebra. The method here used, although perhaps 
less elegant, is, at least, straightforward. 
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with the remainder reflected; etc. Although this infinite sequence of partial absorptions 
and reflections would seem to greatly complicate the analysis, the problem is readily 
solved by the type of method indicated in reference 16. The method will be illustrated 
for the case of the finite -length cylinder that was treated in the preceding section. 

Suppose that the internal surface has an emissivity of and a reflectivity of p i 
(where + p^ = 1). Two simultaneous equations are set up corresponding, respectively, 
to the following two statements: 

(1) The rate u at which radiation leaves a unit area equals the rate at which heat 
enters that unit area through the insulation, plus the rate at which radiation arrives at 
that unit area from the rest of the internal cylinder surface. 

(2) The rate u at which radiation leaves a unit area equals times the local z 
plus p^ times the rate at which internal radiation arrives at that unit area. 

The two statements are translated into the following two equations, respectively: 

(1) u(x) = c[a 0 - z(x)J + i J u( ij) e^ -x d£ + i J u(£) e x "^d£ (22) 

“Z x 

(2) u(x) = € i z(x) + y J u(|) e^' x d£ + y J u(|) e x "^d£ (23) 

If the first equation is multiplied by p i and subtracted from the second equation, the 
integrals cancel and a simple algebraic relation between z and u is obtained: 

u(i - p,)= e t z - Pi c(a 0 - z) 

from which 


z = 


ue i + P iC a 0 
e i + p i c 


Substituting this expression for z into equation (22) gives 


(1 + g)u(x) - a Q g = |J u(£) e^ " x d£ + | J* u(£) e x “^d£ 


(24) 


(25) 


where 


This equation is identical in form to equation (19a). It is solved similarly and has a 
similar solution: 
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(26) 


u = ag - A cosh hx 

where h = \|g/(l + g). The coefficient A is evaluated by the same method as before: 

2a n 

A 2 

e -hZ e hi 

1 + h + 1 - h 

Finally, substituting this derived u into equation (24) gives z, from which T can be 
obtained. 

Some calculations were made for a finite -length cylinder having an internal surface 
emissivity of 0.85. The other conditions were: a = 0.25; e = 0.85; k = 0.0025, 0.025, 
and °°; 1 = 8. The results (not shown) could be compared with the results of the previous 
calculations (fig. 6) for the cases in which all the parameters were the same except that 
e . was 1. As expected, the wall temperatures were found to be slightly higher than those 
for e- = 1. However, the differences were very small, ranging from a few tenths of a 
degree near x = 0 to a maximum of nearly 1° at the open end. Furthermore, the emis- 
sivities of actual mat black surfaces would be considerably higher than 0.85. Accordingly, 
use of unit emissivity was considered to be satisfactory for present purposes. 


Use of the Exact View Factor 


It may be of interest to review a brief study of a different approach to solving the 
simple open tube: The exact view factor (see p. 20) was used in the integral equation, 
but z(x) was represented by a three -term power series, p + qx + rx . Equation (19a) 
is thus replaced by 


(1 + c)( 


p + qx^ + rx 



(x-{)[(x-i) 2 + 6] \ jt i nl 

r 9 ->S/2 ( 2J y 

[(x - zr + iJ 


3/2 


Integrating analytically results in an algebraic equation in p, q, r, and x. Substituting 

three different values of x - say, x = 0, 1/2, and l - gives three simultaneous linear 
equations for p, q, and r, the solution of which gives the power series for z. It is 
here supposed that the three -term series should provide a fairly accurate representation 
of z for a short tube, so that this method should give a fairly accurate result for such a 

tube. 

Some calculations for a short tube (l = 6) showed excellent agreement with the cor- 
responding results in figure 6. For l = 8, the agreement was not quite so close, with a 
maximum discrepancy of 1.5° (at x = 8, with a = 0.25, e = 0.85, and k = 0.0025). For 
still larger l, the disagreement became progressively worse, and the results became 
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physically unreasonable (as, for example, when the temperatures near x = 0 exceeded 
those for the infinite tube). 

The three -term series is obviously inadequate to represent z(x) for the longer 
tubes. Although the representation could be improved by adding terms in x 6 ,x®, 
the amount of algebra would quickly become excessively cumbersome. If a very exact 
solution is needed, it would then be easier to get access to a large computer and use some 
available program for numerically solving such integral equations. 

The fact that good agreement was found for the shorter tubes (l £ 8), where the 
three -term series is reasonably adequate, indicates that the simple exponential view fac- 
tor is basically satisfactory for present purposes. 


Circumferential Temperature Nonuniformity 

For the infinitely long cylinder, with or without heat pipes, figures 3 and 4 showed 
the circumferential temperature variation for several combinations of insulations and 
thermal -control coatings. For a cylinder of finite length, the circumferential variation 
will be greater, mainly because of the reduced internal temperatures. At the reduced 
temperatures, the internal radiation is weaker and less effective in equalizing the internal 
temperatures; accordingly, the temperature difference between the subsolar area and the 
shaded area must increase before equilibrium is established. An additional, related 
reason, especially for the region near the open end, is that the radiation received from 
any longitudinal element of the cylinder is less than it would be if the element were infi- 
nitely long. 


The following reasoning was used in estimating the increased circumferential non- 

OO 

r — v ^ 

uniformity: In equation (13), the term - ) — — cos n0 represents the radiation 


n=0 


received at any Q from the entire infinite cylinder. If the cylinder is of finite length, 
this term would be too large, and should be reduced by some factor s. Furthermore, 
since the longitudinal distributions of T (or of z) along the elements of the cylinder 
would be very similar for all the elements, the factor s should be approximately the 
same for all points around the tube at the same longitudinal location; that is, it may be 
written as s(x). Accordingly, for the finite-length cylinder the term in the equation 

OO 

should be changed to -s(x) ) — cos n 0. The evaluation of s(x) from the local 

^ 4n 2 - 1 

z(x) (the circumferential mean oT^ given, say, by eq. (20)), and its use in estimating 
the circumferential temperature nonuniformity are outlined in the following sections. 
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Without heat pipes.- Equation (13) is modified by multiplying the summation on the 
the left-hand side by s, as just stated. After the term for circumferential heat conduc- 
tion is added, the analysis proceeds precisely as before and yields the following Fourier 
coefficients, which are now functions of s (the subscript f has been used to indicate 
that the series is for a finite -length tube): 



(For s = 1, these coefficients reduce to the coefficients (14) for the infinite cylinder.) 

The coefficient a Q f is still the circumferential mean zj (or oT 4 ), but it varies 
with s. 

As an example of the use of these equations, suppose it is desired to find the cir- 
cumferential temperature distribution at x = 8 radii on a cylinder of length l equal 
to 12 radii, for a ■= 0.25, e = 0.85, and k = 0.0025. From figure 6, the mean tempera- 
ture T at x = 8 is read as 128 K, from which the local mean zi is a x 128 4 = 

3.64 x 10" 4 . Putting this value equal to a 0 f in the first of equations (27) and solving 
for s gives s = 0.981. With this value of s, all the other Fourier coefficients (27) 
may then be evaluated, whence the circumferential distribution of z i and hence of T 
can be determined. 

A number of circumferential distributions were calculated, but they are not pre- 
sented here. Instead, the difference, T max - T min , between the maximum and minimum 
temperatures around the circle, as obtained from these calculated distributions, are 
plotted as functions of T in figure 7 for a = 0.25, e = 0.85, and k = 0.025 and 0.0025. 
It is seen that with decreasing temperature, T max - T m i n increases rapidly. Thus, 
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for the preceding ^example, with T = 128 K, T max - T min is 1.4 K, whereas for the 
infinite cylinder (t = 219 k), T m ax - Tmin is only 0.33 K. For the poorer insulation 
(k = 0.025), the values of Tmax - T m i n become so large at the lower temperatures 
that the basic assumption (the simple modification of eq. (13)) is no longer accurate, 
so that the large derived values of T max - T m i n must be considered as only rough 
approximations. 

With heat pipes . - The preceding analysis can be further modified for the case of 
the cylinders with heat pipes. One first modifies equation (16) by multiplying the second 
and third terms on the left by s. The Fourier coefficients are then determined as 
before: 
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The subsequent treatment is the same as for the case without heat pipes. 

Calculations for a few cases showed that the increase in T max - T • over that 
for the infinite cylinder was relatively small. Whereas for the cal cu late d^ases without 
heat pipes, T max - T min increased by large factors (up to 15) relative to the values for 
the infinite cylinder (fig. 7), calculations for the comparable cases with heat pipes showed 
increases of only about 25 percent. Thus, the isothermalizing effect of conduction to or 
from the heat pipes remains the dominant factor. Calculations for a similarly insulated 
and painted 15 -cm -radius cylinder without heat pipes showed that also for that case the 
isothermalizing effect of conduction restricted the increase to only a few percent. 
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Use of a Heated Isothermal Collar To Provide 
Longitudinally Uniform Temperature 

Basic equations. - The sketch shows again an open-ended tube, except that near the 
end, along a section between x = m and x = n, the wall is maintained at a slightly ele- 
vated uniform temperature by means of heat input from some external source. We sup- 
pose, for the present, that the tube axis is normal to the sun's rays and that the insulation 
is uniform. Along the isothermal section, let oT^ = q; and along the inner and outer 
sections, where the temperatures, in general, are not constant, let oT^ = z(x) and y(x), 
respectively. (See sketch (i).) 


= y(x) 

-l -n -m 0 x m n l 

Sketch (i) 

The determination of z(x) and y(x) involves the solution of two simultaneous 
integral equations. The procedure will be only outlined, since it is not very different 
from what has already been presented. The equation for z(x) is 

z(x) = c[a 0 - z(x)J + | J z(£) e^ _x d£ + | j z(£) e x '^d£ + | qe x ^d£ 



1 r- m 

2 J_ n 


qe«- x d 5 H.I f 
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yU) e^' x d| 


(29) 


When this equation is differentiated twice with respect to x, it yields the differential 
equation 



1 + c v 


- a o) 


Since z is continuous at x = 0 and also must be symmetrical with respect to x = 0, 
it must be of the form 


z = aQ + A cosh fx 



(30) 
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The integral equation for y is 

y(x) = cja 0 - y(x)] + i r" y(« e«- x d? + | y(fl | + § £ y«> e X ' ? d? 


+ | f m qe^" x d£ +i f qe^" X d|+|J z(^) e^' x d£ 

2 J_ n 2 d m * J -m 


(31) 


After differentiating twice, it yields 

SL-X - — £ — (y - ao) 
dx 2 ! + 0) 

Although y is symmetrical with respect to x = 0, it is not defined over the inner region 
(x = -n to n); hence, the solution must be assumed in the more general form 


y = a 0 + Be f|x| + Ce- flx| 


(32) 


In order to solve for A, B, and C, set up three simultaneous equations: 

(1) Substitute expressions (30) and (32) for z and y into equation (29) and per- 
form the indicated integrations for x = 0. 

(2) Similarly, substitute expressions (30) and (32) into equation (31) for y and per- 
form the indicated integrations for x = n and for x = l. 

The three simultaneous equations in A, B, and C that are thus derived are here set 
down since this configuration is of special interest: 

From z(0), 


A<I 


e (f-l)m e -(f-l)m 


f - 1 


f + 1 


+ B 


e (f-Di . e (f-Dn 


+ C 


e -(f+l)n _ e -(f+l)£ 


f-i j [_ 

f + 1 

e - m + e-‘ + e" n ) 

/ -n 

+ q( e - 


From y(n), 


fe-» 

sinh(l + f)m 

\ 2 

1 + f 


1 - f 


+ B r^sh_n[ e (t-l)n . e (i-l«] . (1 + c)e fnj 

+ C j^n[ e -(l+«n _ e -< 1+f »J - (1 + c)e- fn j 

= aQ(e" n sinh n + e'^cosh n - e~ n sinh m) + qe" n (sinh m - sinh n) (34) 
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From y(Z), 


sinh(l + f)m sinh(l - f)m 
1 +f 1 - f 


+ B 


,(f-l)n e (f-l)Z p (l+f)Z _ p (l+f)n 
- + 


e'" - e '‘ ■ / ~ ■ e - (1 + c)e (1+f)i 


2(1 - f) 


2(1 + f) 


+ C 


-(l+f)n _ e -(l+f)l p (l-f)2 _ p (l-f)n 


2(1. + f) 


2(1 - f) 


- (1 + c)e 


a -ml 


= aQ(cosh l - sinh m + sinh n) + q(sinh m - sinh n) 


(35) 


For a given geometry and set of parameters (m, n, l, ag, c (hence f), and q), one 
can evaluate the coefficients of A, B, and C in these three equations and then solve 
simultaneously for the numerical values of A, B, and C. 

Uniform -temperature inner section .- A more interesting problem is to set A = 0 
and solve for B, C, and q from equations (33), (34), and (35); that is, determine the 
temperature (since q = oT 4 ) of the heated collar that will provide a uniform temperature 
between x = 0 and x = m. (The uniformity corresponds to having set A = 0.) For 
this case, since equation (30) simply reduces to z = a q (0 < x < m), the uniform tempera- 
ture is the same as would exist in a similarly painted infinite cylinder. With B, C, 
and q thus known, calculation of the heat that must be put into the collar in order to keep 
it at the required temperature is straightforward. At any longitudinal location on the 
heated collar, the required heat input per unit area is the heat transmission outward 
through the insulation c(q - a Q ] plus the heat radiated inward q minus the heat 
received by radiation from the inside of the cylinder. This last item must be determined 
as the sum of exponential type integrals along the five sections of the cylinder (-1 to -n, 
-n to -m, -m to m, m to n, n to l, where the integral from m to n must be 
in two parts — m to x and x to n). The derived expression for the required rate 
of heat input per unit area must finally be integrated over the collar. 

Calculated results.- A number of calculations were made of the temperature dis- 
tributions and of the corresponding heating powers for a tube of length l - 12 radii with 
a uniform -temperature collar between x = 6 and x = 8 (that is, m = 6 and n = 8). 

The results are shown in figure 8. (Fig. 8(a) is for a - 0.25, e = 0.85, k = 0.0025; and 
fig. 8(b) is for a = 0.15, e = 0.85, k = 0.0025.) The most significant numerical results 
are given in the tables in the figures. (The results for total heat power input to the collar 
are for a tube radius of 150 cm.) 
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The difference between the wall temperatures at x = 0 and at x = 6 may be taken 
as a measure of the temperature uniformity in the 6 -radius -long section between the 
mirror and the collar. For the "design conditions" (the bottom lines of both tables), the 
temperatures at x = 0 and at x = 6 are equal and about 2° below the collar tempera- 
ture. In the worst case (the third line of the table in fig. 8(a)), the telescope was assumed 
to be completely in shadow, or alined with the sunlight, so that there was no external radi- 
ation input. For this case, the temperature at x = 6 is 2.4° higher than that at x = 0, 
and is 3.6° cooler than the collar. 

The 300- to 400 -watt input required to maintain the collar at the specified tempera- 
ture seems relatively modest; and extending the length of the shield (the section from n 
to l) would reduce it to lower values. Furthermore, relatively small amounts of addi- 
tional power applied in the section between the mirror and the collar could maintain a uni- 
form and constant temperature along that section regardless of the tube orientation rela- 
tive to the sunlight. 

The heat source for the collar is a rather challenging design problem. One might 
use some arrangement with variable -conductance heat pipes leading from an isotope power 
unit or from solar-heating panels provided with some phase-change material to furnish 
heat during the shadow phase of the orbit. In any case, some small additional electric 
power would seem to be necessary if very fine control of the temperature is needed. 

Tube With Nonuniform Insulation and Paint 

In all the preceding analyses and examples, the thermal -control coating and insula- 
tion blanket, as characterized by a, e,and k, were assumed to be perfectly uniform 
along the tube. Actually, as previously indicated, the thermal conductance of a super- 
insulation blanket is sensitive to relatively minor irregularities in construction or 
application - for example, the closeness with which the layers are packed. Also, differ- 
ent batches of paint might show some variations, especially in the rate of darkening dur- 
ing exposure to the ultraviolet and particle radiation in space. Some calculations were 
therefore made to determine the effect of changing the values of a and k over a given 
length of an infinite cylinder normal to the solar rays. The simplified geometry was 
considered to be adequate for demonstrating how the effect of an irregularity is washed 
out along a finite cylinder. 

The geometry and definitions are indicated in sketch (j). The subscript 1 is used 
to indicate properties along the variant section between x = -l and x = Z; the subscript 
0 is used to indicate the properties of the remainder of the cylinder that extends to 
infinity in both directions. 
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The mathematical analysis will not be given, since it is similar to the analyses of 
the preceding problems. The solutions are of the form 


Zj = a^ + A cosh f jx 

-fn|X| 

Zo = a Q + Be u 

where 


«0 b 


a ° TC 0 ’ 


ofib 


a l “ W, 


(|x| < l) 
(Ixl > l) 





As before, A and B are determined by solving two simultaneous equations, as for 
z-^O) and z Q (i). 

Figure 9 shows some of the calculated temperature distributions. For all cases, 
ckq was taken as 0.15 and was taken as 0.25. The difference is considerably more 
than should exist in an actual case. The value of e was assumed to be 0.85 for all 
cases, since e does not change appreciably from one batch of paint to another and is 
not much affected by radiation. In some cases, the value of kj was taken as twice the 
value of k 0 ; in others, kj and kQ were equal. A considerable range of l values 
was used. 


For the uninsulated conditions, the curves show large discontinuities, the tempera- 
tures rapidly approaching those for the infinite cylinders at short distances from the dis- 
continuities. For all cases in which there was any insulation at all, the effect of the 
changed a or k is distributed over a considerable length of the cylinder. There is 
a small discontinuity (0.1° to 1°) at the edge of the variant section, which would be some- 
what smoothed out in a practical case, not only because of wall conduction but also 
because changes in insulation characteristics would probably be fairly gradual rather 
than discontinuous. Admittedly, the a and k variations used in these examples are 
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not very realistic. The a variation is probably unreasonably large; and the k varia- 
tion, although not exaggerated, would probably be randomly distributed. Nevertheless, 
it seems clear that one cannot hope to predict temperatures and temperature distributions 
very precisely for a range of conditions simply on the basis of nominal paint and insula- 
tion characteristics or on the basis of model tests. Precision temperature control will 
hardly be achieved without some active controls such as are referred to in the preceding 
section. 


On Maintaining Uniform and Constant Temperature 

Thermal distortions might occur if the telescope temperatures are allowed to 
change, not only because different parts of the telescope may have different thermal - 
expansion coefficients but also because the different parts would change their tempera- 
tures at different rates. It seems advisable, therefore, to strive not only for uniformity 
but also for constancy of temperature. 

It was previously shown that a heated collar near the open end of the telescope tube 
can, under ideal conditions, provide a uniform temperature along the main part of the 
tube (between the mirror and the collar). The uniform temperature in this case is that 
which would be acquired by a similarly painted and insulated infinite cylinder at the same 
orientation relative to the incident radiation. If the orientation of the tube relative to the 
radiation is changed, a different collar temperature could then provide the same mean 
temperature along the tube; however, the temperature will no longer be uniform along the 
tube, although the nonuniformity may not be very large. Thus, a heated collar alone can- 
not suffice to maintain the tube temperature both uniform and constant as the orientation 
(hence the radiation input) is varied. 

One approach to uniformity and constancy might be to use the highest temperature 
condition (the tube axis normal to the solar rays) to define the design temperature, and 
add additional heat along the tube to maintain this uniform temperature for all other 
orientations. For k = 0.0025, this additional heating for the 150-cm-radius tube would 
be only about 50 watts for the worst case (no external radiation input at all, as with the 
tube axis parallel to the solar rays) . 

Longitudinal heat pipes .- A method of maintaining a more uniform temperature 
along the tube would be to provide a set of longitudinal heat pipes. If, in addition, these 
longitudinal heat pipes were thermally connected together, they would also contribute to 
improving the circumferential temperature uniformity. In reference 4 it was suggested 
that a shallow drumlike vapor chamber might be located just under the mirror, and that 
the longitudinal heat pipes might be brought into the vapor chamber, as indicated in 
sketch (k). Such an arrangement would provide not only a uniform radiation environment 
for the bottom of the mirror, but would also insure that the temperature of this radia- 
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tion was very nearly the same as that of the telescope tube itself. However, actual con- 
struction of a flightworthy vapor chamber of this type and size with many inserted heat 
pipes seems a rather formidable undertaking. Perhaps merely bonding the heat pipes to 
the outer top or bottom surface of the vapor chamber would be adequate. 

Balloon-type enclosures.- In reference 1 it was suggested that the thermal problem 
could be eased by enclosing the entire telescope in a balloon. The thermal -radiation 
environment within such a balloon should be nearly uniform and isotropic and should also 
be nearly independent of orientation relative to the external radiation, provided the inner 
surface of the balloon is nearly specular. Reference 17 contains some space -simulation - 
chamber studies of the radiation inside a model balloon irradiated by a solar simulator. 
On the whole, the data tend to confirm the prediction. For example, in one case, the 
external surface had a mean a/e ratio of 1.15 (for which the mean 
internal temperature was 290 K), the internal surface was lightly 
anodized, and the outer surface of the "telescope tube" was specular. 

The temperature of the radiation striking the outside surface of the 
telescope tube (sketch (i)) was found to be uniform to within about 
±5 K for a range of orientations. This uniformity, considered in 
terms of (radiation temperature)^, is over 20 times better than for 
an unenclosed tube irradiated from one side by sunlight. Furthermore, the mean temper- 
ature of the radiation incident on the tube was, as expected, almost independent of the 
orientation of the tube relative to the solar -simulator beam. 



Thus, the balloon could be very helpful in providing relatively constant and circum- 
ferentially uniform temperatures. Furthermore, if the balloon were extra large, with a 
diameter considerably exceeding the telescope length (sketch (m)), the longitudinal non- 
uniformity, with decreasing temperatures toward the open end (as 
shown in fig. 6), could also be greatly reduced. (Such a spacious 
volume around the telescope was probably intended in ref. 1; in 
fact, it mentioned a set of concentric balloons.) More extreme uni- 
formity and constancy would require additional complexity, such as 

the previously discussed vapor chambers, heat pipes, heated collars, Sketch (m) 
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fine -control heating power, and insulation; however, the performance requirements of all 
such components would now be greatly reduced. 

Thus, a well -insulated telescope enclosed within an extra-large balloon and provided 
with perhaps 150 watts (for the 150-cm-radius telescope) of distributed and regulatable 
heating power appears to be an attractive concept. Certain inherent elements of awkward- 
ness in the concept have probably prevented its serious consideration in such systems 
studies as that of reference 18. For example, the balloon could not be a simple sheet, 
like the Echo satellites, but would have to be sufficiently stiff, probably with internal 
bracing, so that it could be readily reoriented as required. The large solar-cell panels 
extended from the side of the telescope tube in the design proposed in reference 18 would 
have to be located outside of the sphere, somehow attached to the bracing. A torn area, 
such as might be produced during deployment, might not seriously affect the thermal con- 
trol, but might let in sunlight that would be reflected around the inside of the balloon and 
into the telescope. During the shadow phase, the temperature inside the balloon would 
drop, so that consideration would still have to be given to means of preventing periodic 
variations of telescope temperature in response to the alternating sunlight and shadow of 
a low orbit. 


LOCAL WALL HEATING VARIATION DURING ORBIT 

If the telescope is in a fairly close orbit about the earth, it will experience alter- 
nating phases of sunlight and shadow, with a period of about 100 minutes; or, if it is in a 
geosynchronous orbit, it may experience daily shadow periods lasting as long as 1 hour. 
The question now arises as to whether a fairly well insulated telescope would experience 
any significant heat -input oscillations in response to the unsteadiness of the thermal input 
to the outer skin. 

Any analysis of this problem for a relatively undefined structure, and without a 
thorough knowledge of the characteristics of unsteady heat flow through multilayer insu- 
lation would seem to carry a rather heavy burden of assumptions. Some experiments in 
a space -simulation chamber will be needed in order to procure the basic quantitative 
data. For present purposes it was considered that a simplified analysis would suffice to 
provide some insight into the problem. 

Consider that the insulation consists of 50 aluminized plastic sheets, each 
0.00075 cm thick. A "radiation conductance" of 0.0025 was assumed, as in many of the 
calculated examples. However, this value corresponds theoretically to only about 
10 sheets, since the radiation conductance between a pair of adjacent sheets is about half 
the emissivity of the aluminized plastic sheet, or 0.025. Accordingly, the 50 sheets were 
represented, for the analysis, by 10 sheets, each five times as thick as an actual single 
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sheet, or 0.00375 cm thick; and the radiation conductance between adjacent pairs of sheets 
was retained as 0.025. 


The temperatures of the 10 sheets, starting with the outside sheet, are designated 
T 0 , Tj, . . T 9 ; and the temperature of the cylinder wall is designated T 19 . The wall 

temperature T^q is assumed to be controlled to a constant value, equal to that for an 
infinite tube with the given a and e oriented normal to sunlight. The outside sheet 
quickly attains radiation equilibrium with the environment, so that its temperature is 
assumed to be known as a function of time. One may then set up nine equations for the 
other temperatures Tj, T2, . . Tg: 


0.025a (t q 4 
0.025a(T 1 4 


+ T 2 4 - 
+ T 3 4 - 


2Tj 4 ) x 0.00375 X 0.4 
2T 2 4 ) = 0.00375 X 0.4 


dT x 

"dT 

dT 2 


(34) 


0.025a(T 8 4 + T 10 4 - 


2T 9 4 ) = 0.00375 X 0.4 


dt 


where the plastic is assumed to have unit density and a specific heat of 0.4. 

Calculations were made for four cases. For each case, the nine temperatures 
Tj, T 2 , . . ., T 9 were determined as functions of time. However, only T9 (the tem- 
perature of the sheet immediately adjacent to the wall) will be shown, since it deter- 
mines directly the rate of radiation heat input to the wall from the outside. The calcula- 
tions were made only for a subsolar area. 


Synchronous Orbit 

For a synchronous orbit, the subsolar area was assumed to have an outside tem- 
perature Tq that is constant at 396 K except for 1 hour when it is essentially zero. 

The wall temperature Tjq was assumed to be constant at 297.4 K (corresponding to 
a = 0.25, € = 0.85). Figure 10 shows the response of T9 to the 1-hour shadow phase. 

Two curves are shown. The curve with the deeper dip is for 50 layers of aluminized 
plastic; the other is for 100 layers. The analysis for the latter case merely required 
that the factor 0.00375 in equations (34) be changed to 0.0075 (since each of the 10 sheets 
is now considered to be ten times as thick as an actual single sheet). 

The lag and the spreading out of the response are very apparent. For the thinner 
insulation, there is essentially no response for the first 15 minutes after the beginning 
of the shadow phase, the minimum is reached about 1/2 hour after the end of the shadow 
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phase, and the entire perturbation lasts about 4 hours. For the thicker insulation, the 
initial lag is about 30 minutes, the minimum is less deep and occurs about 1 hour after 
the end of the shadow phase, and the entire perturbation lasts about 7 hours. 

Spreading the response over several hours is a useful function of the heavy multi- 
layered insulation. The solar heat input to the wall of the entire telescope tube may be, 
say 40 watts; so during the 1 -hour shadow phase, a 40-watt -hour deficit occurs, which 
must be supplied from other sources - probably electrical - if the wall temperature 
is to remain constant. The fact that the response is spread over, say, 7 hours means 
that an average power of only 6 watts is needed during this period. For exact compensa- 
tion, however, about 20 watts would be needed for the period when Tg is near its mini- 
mum (or 32 watts for the thinner insulation). 

The time between emergence from shadow and return of Tg to its normal value 
(about 3 hours for the 50-layer blanket and about 6 hours for the 100-layer blanket) may 
be taken as a rough indication of the time required after a reorientation of the telescope 
for Tg to attain its new equilibrium value. 

Low Orbit 

For the low orbit, alternating sunlight and shadow phases of 60- and 40-minute 
durations, respectively, were assumed. The wall temperature T 10 was again taken as 
297.4 K, and, for the subsolar area under consideration, the temperature of the outermost 
layer T 0 was assumed to be 396 K when in sunlight and 244 K when in shadow. This 
variation, shown in the upper curve of figure 11 , is obviously a gross simplification of 
any actual case, but it should be adequate for present purposes. The calculations were 
made, as for the preceding problem, for 50-layer and 100-layer insulation blankets. 

The lower curves of figure 11 show the calculated variation of Tg with time. 

The oscillation of Tg lags that of T 0 by about 45 minutes for the 50-layer blanket, 
and by about 75 minutes for the 100-layer blanket. The amplitude of the oscillation is 
four times as much for the 50-layer blanket as for the 100-layer blanket. Thus, the 
heavier blanket provides a strong damping of the oscillating input, but the oscillation of 
T 9 for the lighter blanket is a fairly large fraction of the mean difference between Tg 
and the wall temperature Tio- Furthermore, since the amplitude seems to vary 
inversely as the square of the number of layers, it would probably increase rapidly with 
further decrease in the number of layers. 

RESUME AND CONCLUDING REMARKS 


Analytical methods have been developed for estimating the temperature distribution 
in a laige space telescope, idealized as a simple insulated cylindrical tube, open at one 



end, and with a flat perfect mirror across 1 the other end. Since heat transfer, both 
within the tube and through the insulation, is mainly by radiation, oT^ was taken as the 
basic variable (where <j is the Stefan-Boltzmann constant and , T is the temperature). 
Expressing the circumferential distribution pf °T^ as a Fourier series permitted a, 
straightforward determination pf the circumferential distribution; and use of an exponen- 
tial expression to approximate the view factor between circular, bands of the tube pro- 
vided tractable equations for the longitudinal temperature distribution. Methods for 
taking into account the effects of multilayer insulation, thermal-coating characteristics, 
heat pipes, and heated collars were also developed. • . ... ... . 

A few numerical calculations were made to show The ^significance- of .each of these 
components. The cylinder radius was taken as 150 cm for most of the calculations. 

Where wall conduction was involved, the wall was assumed to be of 0.1 -cm -thick alumi- 
num. Where circumferential heat pipes were considered, they were assumed to be 40 cm 
apart. The cylinder axis was taken as normal to the solar rays, and the shadowed half 
was assumed to receive no radiation input at all. These assumptions can be altered with- 
out significant modification of the analytical methods. 

If the tube is wrapped in a multilayer insulation, the circumferential temperature 
nonuniformity can be limited to tenths of a degree. If, in addition, circumferential heat 
pipes are attached, the nonuniformity can be further reduced to values of the order of 
hundredths of a degree. A multilayer insulation also tends to minimize any temperature 
irregularities of the tube wall due to nonuniformity of the thermal -paint characteristics 
or of the insulation itself. It can also damp out time variations due to alternating sun- 
light and shadow of a low orbit. On the other hand, thermal equilibrium, at which the 
heat transmitted in through the insulation equals the heat radiated out of the open end of 
the tube, may exist at undesirably low internal temperatures when such insulation is used; 
also, there will be a large temperature drop along the tube from the mirror to the open 
end. A heated constant -temperature collar near the open end of the cylinder will largely 
eliminate both of these problems. 

Although this paper is largely concerned with temperature uniformity of the tele- 
scope, just how much uniformity is essential or desirable remains an open question. The 
development of glasses (and perhaps also structural materials) with very low thermal- 
expansion coefficients, and the development of means for adjusting the mirror figure 
(known as "active optics”) would seem to ease the thermal -control requirements. The 
availability of these materials and techniques, however, would hardly justify designing 
for minimal thermal control. For example, even if active optics is used, it would never- 
theless be preferable if checking and adjusting the optics were required only occasionally, 
rather than every few minutes. Other types of technical developments might also appre- 
ciably influence the requirements. For example, the main structure between the primary 
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tind secondary mirrors might be a tubular framework, constructed as one continuous 
heat pipe. 

In any case, the purpose of the present paper is not to propose or advocate any 
particular degree of temperature control or any particular means of attaining it, but 
rather to present some practical analytical methods and some calculated results for a 
first -approximation configuration, 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 16, 1973. 
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Figure 2.- Measured conductivities of several multilayer insulations, from references 9 
and 10. Best -fit straight lines are shown for every set of data. The heavy solid line 
has a slope of 3 and corresponds to a "radiation conductance" k of 0.0025 for a 
1-cm-thick layer (conductivity, 4oT^k). See line key on page 45. 
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KEY FOR FIGURE 2 


Symbol 

A 

□ 

o 

o 

0 

A 

O 


Construction Packing density 

(layers and separators) (layers/cm) 

Alternate layers of dimpled singly 75.5 

aluminized Mylar and embossed 
singly aluminized Mylar 

Doubly aluminized Mylar 9.1 

(perforated) - red foam 

Doubly aluminized Mylar - red 9.1 

foam 


Embossed singly aluminized 19.7 

Mylar - nylon net 

Doubly aluminized Mylar- 8.7 

fiber-glass tufts 

Crinkled singly aluminized Mylar 27.6 

(NRC-2) 

Embossed singly aluminized Mylar 63.8 
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Figure 3.- Temperature distribution around infinite cylinder with axis normal to solar rays. Wall conduction taken 
into account. Cylinder wall, 0.1-cm-thick aluminum; radius, r = 150 cm unless otherwise noted. 
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(a) No insulation. 
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(b) Insulated, k = 0.0025. 


Angle from subsolar point, deg 

Figure 4.- Temperatures around infinite cylinder with heat pipes at 40 cm. Axis normal 
to solar rays; radius, 150 cm; wall, 0.1-cm-thick aluminum. Temperatures are on 
the circle midway between adjacent heat pipes. 
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Figure 6.- Variation of temperature (circumferential average) along inside wall of 
insulated and uninsulated open-ended tubes. Tube axis normal to solar rays. 
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(e) a = 0.25; e=0.85 ; k = 0.0025. 


(g) a = e = 0.4; no insulation. 
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Figure 6.- Continued. 
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Figure 6.- Concluded. 
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Figure 7.- Difference between maximum and minimum wall temperatures around the 


cylinder as a function of local root-mean-fourth-power temperature T 
Finite-length tube with axis normal to solar rays; a = 0.25; e = 0.85; 


(=(a 0 /a) 1/4 ). 
k = 0.025 


and 0.0025; r = 150 cm; no heat pipes. 
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(b) a = 0.15; e = 0.85. 
Figure 8.- Concluded. 
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Figure 9.- Temperatures along wall of infinite cylinder having a short section with altered 
properties of the insulation, paint, or both. Cylinder axis normal to the solar rays. 
The short lines at the right indicate the asymptotic temperature - that is, the temper- 
ature in an infinite cylinder with the given values of ap and e Q . The short lines at 
the left indicate the temperature in an infinite cylinder with the given values of a j 
and e i . 
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Figure 10. Response to a 1-hour shadow phase of temperature of sheet adjacent to tube wall. Subsolar area 
tube axis normal to solar rays; a/e = 1; k = 0.0025; wall temperature, 297.4 K. 
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Figure 11.- Response to alternating 60-minute sunlight and 40-minute shadow phases of 
the temperatures of the outermost sheet (To, assumed) and of the sheet adjacent to 
the tube wall (Tg). Wall temperature, 297.4 K; k = 0.0025. 
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